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Velocity oscillations in confined channel flows of concentrated colloidal suspensions 
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We study the pressure-driven flow of concentrated colloids confined in glass micro- channels at the 
single particle level using fast confocal microscopy. For channel to particle size ratios a/D < 30, 
the flow rate of the suspended particles shows fluctuations. These turn into regular oscillations for 
higher confinements (a/D ~ 20). We present evidence to link these oscillations with the relative 
flow of solvent and particles (permeation) and the effect of confinement on shear thickening. 

PACS numbers: 83.80.Hj, 83.50.-v, 83.60.Rs, 83.60.Wc 



Confined flows are important in both fundamental 
and applied science. In simple liquids, confinement ef- 
fects show up on the ten- molecule, or ~ 10 nm, length- 
scale p], Q. The properties of soft and granular mat- 
ter, however, are controlled by mesoscopic objects (poly- 
mers, colloids, grains), whose sizes can be fim to mm 
or beyond. Here, confinement effects can be relevant on 
length scales that are macroscopic, e.g. in the silo flow 
of grains [3], or at least on length scales resolvable in 
the optical microscope, such as the flow of suspensions 
of micron-sized particles in ~ 10 — 100/im devices. Thus, 
for example, confinement effects in the nozzle-directed 
printing of complex structures using colloidal inks [4|, [5[ 
have recently been found p . More generally, colloids are 
increasingly found in microfluidic applications [7j. 

Here we present experiments on the flow of a hard- 
sphere suspension at nearly random close packing, a 
'paste' confined in 20 to 40 particle diameter square 
channels. Strong confinement is found to affect the qui- 
escent properties of dense colloids, e.g. the onset of the 
hard-spheres glass transition [9], and therefore we may 
expect strong effects under flow, too. By tracking parti- 
cles using fast confocal microscopy, we indeed find large- 
amplitude velocity oscillations in the narrowest channels. 

Oscillations are ubiquitous in paste flow [TO, EH, G3 , 
but their origin remains unclear. The oscillations we see 
are clearly due to confinement — they are prominent in 
20-diameter channels, but vanish for 40-diameter chan- 
nels. Interestingly, some channels used in some previous 
experiments, when expressed in particle diameters, are of 
the same order of magnitude as those used in our work: 
~ 10 — 40 particle diameters in [l(| and ~ 50 — 100 di- 
ameters in [ll|; it is possible that confinement may also 
be relevant in these cases. 

We used sterically-stabilized polymethylmethacrylate 
(PMMA) spheres (diameter D = 2.6 ± 0.1 /xm; vol- 
ume fraction $ > 0.63, from confocal microscopy; Brow- 
nian time tb = 4.4s), fluorescently labelled with ni- 
trobenzoxadiazole and suspended in a mixture of cyclo- 
heptylbromide and mixed decalin (viscosity 2.6 mPa-s) 
for buoyancy matching at room temperature. A pres- 
sure difference, AP (« 10 2 to 10 4 Pa), is applied to 
drive the suspension into square borosilicate glass micro- 
channels (Vitrocom Ltd; width a = 50, 80 and 100/im) 
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FIG. 1: (a) Normalized particle volume fraction as a func- 
tion of x from particle counting in 1000 frames, (b) Velocity 
profiles as a function of normalized y (a e ff = a — D/2) cor- 
responding to points 1 and 2 in (a), (c) Velocity offset by 
the long-time average as a function of time for three different 
channel sizes, (d) Space-time (x-t) diagram of oscillating flow 
in a ~ 20D channel. The arrows highlight jamming events. 



with untreated, smooth inner walls. The channel dimen- 
sions (a/D ~ 20, 30 and 40) were not integral multiples 
of D to avoid complete layering. 

The flow across the full width of the channel y was 
imaged with a Visitech VTeye confocal scanner and 
a Nikkon TE 300 inverted microscope. We collected 
44/imx58/im images (107 frames/s) at 17 /urn from the 
lower surface. From these we located (resolution ~ 
50 nm) and tracked particles [l3|, 0]. From the sin- 
gle particle velocities we calculated the macroscopic flow 
velocity, averaged over five frames to eliminate high fre- 
quency instrumental and tracking noise. 

In the widest channel (a/D ~ 40), the velocity is 
steady with time, Fig. [jjc [24]. The velocity profile and 
density (data not shown) do not vary with position along 
the channel x after entrance effects. The velocity pro- 
file consists of shear zones near the walls and a central 
plug with a flow-speed-independent size. We have re- 
cently shown that such profiles can be understood using 
a model developed for granular pipe flow [151 ]. 



2 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

f[Hz] 

FIG. 2: Normalized power spectrum for the velocity signal 
in the inset to Fig. [^((y) = 20.4/im/s). Inset: four seconds 
of the offset velocity signal with the fit obtained from the 
superposition of 6 sine waves from the power spectrum (red). 
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FIG. 3: Flow velocity V(t) (black, left axis), offset by the 
long-time average (V). Right axis: corresponding mean (V) 
relative to the threshold for oscillations (normalized by Vthr = 
7/xm/s, red bullets), obtained by changing the applied AP 
in the highlighted time windows (veritcal segments). Inset: 
example velocity oscillations, with one cycle highlighted (red). 



The situation is radically transformed when the chan- 
nel size is halved (to a/D ~ 20). The average veloc- 
ity profile and density are now x-dependent, Fig. [l]-a. 
Near the inlet (x = 0) we find the same kind of sheared 
profile as in the wider channels. Further downstream 
we encounter another region where the central plug ex- 
tends across the full channel, i.e. shear between parti- 
cles is completely absent and the suspension flows as a 
solid body with velocity V; in this complete plug region 
the density reaches a maximum value <£ p , and remains 
constant over ~ 1cm, before terminating even further 
downstream at a sharp interface with another sheared 
region, Fig. [T]-a, whose density decreases away from the 
plug. Dramatic, quasi-regular velocity oscillations are 
also seen, Fig. [l]-c. Such oscillations are nearly washed 
out in a ^ 30D channel, and absent when the channel is 
widened to ~ 40D. 

Fig. [l]-d shows a 'space-time diagram' of the oscil- 
lations in the c± 20D channel, constructed by stack- 
ing one-pixel-wide ?/-bins taken from consecutive images 
in the channel's center, in the region of complete plug 
flow. Changes in the slope of the traces correspond to 
changes in the flow speed, and we observe the presence 
of regularly spaced horizontal bands corresponding to 
events where the flow comes to an almost complete arrest 
(jams). The power spectrum of the oscillations, Fig. O 
consists of a fundamental frequency and higher harmon- 
ics; superposing the fundamental and the first five har- 
monics reproduces the gross features of the signal (inset). 

These regular jams are absent below a particular 
threshold average flow rate Vthr- Fig.[3]shows V— (V) and 
the corresponding normalized distance from the thresh- 
old [(V) — Vthr]/Vthr versus time (a ~ 20D). For small 
(V) the flow is time-independent, but above the thresh- 
old oscillations appear. Reducing the average flow rate 
reduces the amplitude of the velocity fluctuations. 



We measured the fundamental frequency of the os- 
cillations as a function of the local average speed (V) 
(Fig. 2]) [25|]. The flow is steady for average velocities 
below a run-dependent threshold, while oscillations arise 
above this threshold (bottom inset, Fig. 2]). The frequen- 
cies / can be fitted by f((V)) = a((V) - V thr f, where 
/ extrapolates to at Vthr (for the curve in the inset 
Vthr = 10.06/im/s) and a and (3 are fitting parameters. 
Plotting / versus (V) — Vthr collapses all the data onto 
one single curve, which is in turns fitted by a power law 
with exponent (3 = 0.34 ± 0.13. 

The amplitude of the oscillations (estimated by the 
standard deviation of the velocity signal, Fig. [4j top in- 
set) shows a linear increase with flow speed above the 
threshold. The degree to which the linear fit misses the 
origin represents the upper limit of the noise in non- 
oscillating signals. This linear dependence is directly re- 
lated to the fact that at each jamming event, the velocity 
drops from roughly « (V) to a much lower value (see the 
'dips' in the inset to Fig. [3]) that is only weakly dependent 
on the average flow speed. 

Summarizing our observations, we find that under 
strong confinement (a/D ~ 20) and above a threshold 
flow rate, the flow velocity oscillates while driven by a 
constant pressure difference. The frequency and the am- 
plitude of the oscillations increase with flow rate. The ve- 
locity profile and particle density are non-uniform along 
the channel. The velocity oscillations and spatial non- 
uniformities along the flow directions disappear when the 
channel width is doubled (a/D ~ 40). 

To understand the physical origins of these oscillations, 
we analyse a typical cycle of the velocity of the plug re- 
gion in detail (black curve, Fig. [SJ. For convenience, we 
start at a time when the plug is jammed (due to a mech- 
anism that we will shortly propose), Point 1 in Fig. [5j 
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FIG. 4: Fundamental frequency of the velocity oscillations as 
a function of (V) — VthrJ the full curve is a power- law fit with 
exponent 0.34 ± 0.13. Top inset: amplitude (standard devi- 
ation) of the velocity oscillations as a function of (V) — Vthr 
(points) and a linear fit (line). Bottom inset: fundamental 
oscillations frequency as a function of the average flow speed 
(V) for one run (points) with a power- law fit (curve). 



The (constant) pressure gradient along the channel will 
drive permeation flow of the solvent relative to the plug. 
We suggest that this erodes particles from the front of 
the plug, and gives rise to a less dense region; such 
a region is indeed observed, Fig. [T]-a. A density wave 
should therefore propagate backwards, slightly expand- 
ing the plug, and allowing the whole plug to accelerate 
from rest, Point 2, Fig. [5l We see evidence for such 
a back-propagating dilatorial wave by analysing parti- 
cle displacements in the plug's co-moving frame. The 
blue curve in Fig. [5] reports the absolute value of the x- 
component of the particle velocity in this frame, |z/ x |. We 
observe a rise as the plug unjams. The (signed) values of 
this velocity are represented in the left-hand inset image, 
and clearly show an expansion of the plug. 

The less dense region at the front of the plug can sus- 
tain shear (point 2, Fig.[I]-a). We estimate that the shear 
rates here (in terms of the Peclet number, Pe) are in the 
region 2 < Pe < 20. In this regime, concentrated suspen- 
sions are expected to shear thicken [16]: a consequence of 
the inability of particles to move past each other quickly 
at high particle densities. We suggest that as the shear 
front begins to thicken and slows down, particles accumu- 
late at its back (upstream), increasing the local density 
and enhancing the tendency to shear thicken, until the 
suspension suddenly jams. 

In this scenario, we expect a compression wave passing 
through the system at this point. The blue curve in Fig. [5] 
does indeed show a sharp peak at Point 3. The (signed) 
values of this velocity represented in the right-hand inset 
also show this compression. The cycle now repeats. 

The existence of a threshold flow rate (Pe c± 2) for 
oscillations can be traced back to a critical rate needed 
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FIG. 5: Plug velocity as a function of time in the lab frame 
(right axis) and absolute value of the x-component of particle 
velocities in the comoving reference frame. The values of z/ x 
have been calculated using particle displacements over 0.1 s. 
The inset diagrams show the spatial distribution of zy x over 
the field of view, with flow direction being from right to left. 



for the onset of shear thickening. Moreover as recently 
found by Fall et al. [17[, confinement plays a crucial role 
in the occurrence of shear thickening, with the critical 
shear rate for shear thickening decreasing as the gap size 
decreases in the 10-100 particle diameter regime. Appar- 
ently, then, within the range of flow rates we can explore, 
the channel needs to be smaller than c± 30Z) for shear 
thickening, and therefore oscillations, to occur. 

The interplay between permeation and shear thicken- 
ing can therefore give rise to oscillations in narrow chan- 
nels. However, the regularity of the oscillations observed 
under strong confinement requires further discussion. 

We have already suggested that permeation unjams the 
system by eroding particles at the plug-shear interface. 
Given that particles close to the centerline move faster in 
the shear region, (Fig. [TJ-b), the streamlines in the plug- 
to-shear transition zone must be convergent (Fig. [(J left), 
reminiscent of particle movement in the vicinity of an im- 
posed geometrical constriction, such as an hourglass or 
funnel. There is, however, an important difference. The 
colloids in our case self- organize into a 'constriction with 
convergent flow' scenario, i.e. a fraction of the colloids 
consitute the constriction into which the rest of the par- 
ticles flows convergently. We suggest that it is such self- 
organization, with appropriate feedback, that creates the 
regular oscillations we see under high confinement. 

If our surmise is correct, then we expect the oscilla- 
tions to become much less regular if the geometric con- 
striction is imposed rather than self-organized. This is 
indeed seen if the channel walls are made rough on the 
particle scale (by sintering a disordered layer of colloids 
onto the channel walls [l5[ ) . Boundary roughness hinders 
the presence of a complete plug and imposes jamming 
at specific ('quenched') locations in the channel. Fig. [6] 
(right) shows a space-time diagram of a flow in a ^ 20Z) 
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FIG. 6: Left: superposition of particles positions during flow 
(right to left, (V) ~ 25/xm/s, 50 images) at the plug-shear 
interface showing converging streamlines. Right: space-time 
diagram of erratic flow in a ~ 20D channel with rough walls. 
The arrows highlight jamming events. 

channel with rough walls. Under an intial applied AP 
the flow jams completly at time to (not in the figure) 
and stays jammed until AP, and thus permeation, is in- 
creased tenfold. The flow then starts again with a series 
of erratic jamming events. Erratic jamming was also re- 
ported by Haw in the proximity of a fixed constrction in 
the flow of sterically-stabilised PMMA colloids [Hj]. 

The mechanism described above is suggestive of sim- 
ilarities with the case of dry grains flowing in an hour- 
glass [18], where the regularity is due to the interplay of 
air counterflow and jamming of the particles at the con- 
striction pjl [2(j. Here, the regularity is inherent to the 
particle dynamics under confined shear flow. The 'self- 



organized' constriction and jamming in our channels ap- 
pear only if the shear rates are sufficiently high, i.e. if 
the suspension is shear thickening. Moreover, for rates 
above the shear thickening threshold, we measured that 
it takes a certain accumulated strain (~ 2 to 10) for the 
system to jam, suggesting analogies with strain harden- 
ing in colloidal suspensions [2l| . The strong confinement 
effects we observe also recall analogies in dry granular 
flow |3, Q. 

In conclusion, we have observed strong oscillations 
in the flow of concentrated hard-sphere-like colloids in 
square channels. Oscillations are widespread in the bulk 
flow of pastes. In previous studies using polydisperse 
industrial particles, oscillations were ascribed to filtra- 
tion [10] and non-specific 'stick-slip' instabilities [ill ]. 
By studying a well-characterised model suspension using 
real-time, single-particle imaging, we have been able to 
propose a specific, microscopic mechanism for the onset 
of instability: an interplay between solvent permeation 
and shear thickening in the plug-to-shear transition zone. 
We observe such instabilities only in channels narrower 
than ~ 30Z), likely because of confinement effects on 
shear thickening [17[. Interestingly, confinement effects 
are probably also relevant in previous studies (~ 10— 40Z) 
[l0[ and - 50 - lOOZ) [n|, but note the much higher 
Pe because of particle size in both cases). Irrespective 
of whether our proposed mechanism is also relevant for 
these industrial particulates, our results should have di- 
rect implications for the flow of concentrated colloids in 
microfluidic devices @, 0] • 
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2-3 times larger than the typical oscillation periods. 
[25] This can also be done in terms of AP, but the former 
represents better the microscopies of the flow, since the 
pressure difference is only controlled macroscopically and 



we are not able to access its local values. Moreover the 
relation between applied AP and measured (V) can be 
non-monotonic and strongly history dependent. 



